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iT u and S* production in Au+Au collisions at ^s NN = 200 GeV and 62.4 GeV 
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Applying a quark combination model for the hadronization of Quark Gluon Plasma (QGP) and A Relativistic 
Transport (ART) model for the subsequent hadronic rescattering process, we investigate the production of K*° 
and E* resonances in central Au+Au collisions at -\/s NN = 200 GeV and 62.4 GeV. The initial K*° produced via 
hadronization is higher than the experimental data in the low p T region and is close to the data at 2-3 GeV/c. We 
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take into account the hadronic rescattering effects which lead to a strong suppression of K with low p T , and 



find that the p T spectrum of A** can be well described. According to the suppressed magnitude of A"* yield, 
the time span of hadronic rescattering stage is estimated to be about 13 fm/c at 200 GeV and 5 fm/c at 62.4 
GeV. The pr spectrum of Z* directly obtained by quark combination hadronization in central Au+Au collisions 
at 200 GeV is in well agreement with the experimental data, which shows a weak hadronic rescattering effects. 
The elliptic flow v2 of E* in minimum bias Au+Au collisions at 200 GeV and p T spectrum of E* at lower 62.4 
GeV are predicted. 
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I. INTRODUCTION 

The short-lived resonances are efficient tools of probing the properties of the hot and dense medium produced in 
relativistic heavy ion collisions. At RHIC energies, QGP with extremely high energy density is created in primordial 
collision stage [ 1], and the system undergoes a long time to expand and cool [2]. The lifetime of the resonance is about 
a few fm/c, which is less than (or roughly the order of) the lifetime of the system formed in heavy ion collisions. After 
QGP hadronization, but before the interactions of hadrons cease, the initially produced resonances and stable hadrons 
will undergo a hadronic rescattering stage. The resonance might be destructed by rescattering with other hadrons and 
also be regenerated by the collisions of other hadrons, and decay daughter particles of resonance are kicked by other 
hadrons causing the signal loss. The physical properties of resonances, e.g. their masses and widths, might be modified 



'Electronic address: shaofl@mail.sdu.edu.cn 



by the surrounding medium |3|,|4j]. In addition, the yields and momentum spectra of resonances might be changed. The 
experimentally reconstructed resonances are the synthetic results of the hadronization and hadronic rescattering effects. 

RHIC and SPS experiments have provided rich data of production of resonances such as K*, X* in relativistic heavy 
ion collisions 15|48[]. The relative yield ratios of resonances to stable hadrons are studied experimentally, which also 



invokes many theoretical explanations ]9 



1111 . This progress greatly promote our understanding of QGP hadronization 



mechanism tested against stable hadrons and the hadronic rescattering dynamics, e.g. the time span of hadronic stage and 



cross sections of various hadronic interaction channels [J12 , 



13] 



The production mechanism of resonances at QGP hadronization is hard to identify due to the entanglement of 
hadronization and rescattering effects. RHIC data, e.g. the phenomena of v2 quark number scaling and high p/n ratio etc, 
show that the production of various stable hadrons at QGP hadronization is realized by the combination of constituent 



(anti-)quarks Ill4l - ll6ll . Recently, STAR experiments observed that the v2 of K*°, the same as stable particles, follows 
the constituent quark number scaling rule in the intermediate pj region [5Q. This provides an evidence for the quark 
(re-)combination/coalescence mechanism of resonance production at hadronization in relativistic heavy ion collisions at 
RHIC. The K* meson andE* baryon are of particular interests due to their very short lifetime ( 4fm/c) and strange valence 
quark content. The experimental data of their midrapidity yields and pj spectra are all available recently |]5l |6J . A system- 
atical study of this pair of meson and baryon resonance should be able to further test the quark combination mechanism 
of the hadron production in relativistic heavy ion collisions. 

In this paper, we apply a quark (re-)combination model for the hadronization of hot and dense quark matter and ART 



model 1117 , 



1811 for the hadronic rescattering process to study the K*° and 2* production in central Au+Au collisions 
at Y^aw =200 GeV and 62.4 GeV. The investigation strategy is divided into two steps. Firstly, we compare directly the 
hadronization results with experimental data to investigate the proportion/magnitude of hadronization exhibited in the final 



20J] to 



observation. Here, as a tool, we use the quark combination model developed by Shandong group (SDQCM) 11191 
deal with QGP hadronization. Secondly, we take into account effects of hadronic rescattering and study the entanglement 
of hadronization and hadronic rescattering effects at RHIC 200 GeV and 62.4 GeV and time span of hadronic stage for 
the system produced at high RHIC energies. 



II. INITIAL K'° PRODUCTION VIA HADRONIZATION 



The performance of quark (re-)combination mechanism on explaining the 



intermediate pj region in relativistic heavy ion collisions is pretty well IU4M16 



>roduction of various stable hadrons in the 
2111 . The mechanism can also well describe 



Pt -integrated yields and rapidity spectra of hadrons at RHIC and high SPS energies [2C , 



recombination models at RHIC. Quark recombination model 



describe the combination of quarks into hadrons. ALCOR I125I1 and SDQCM apply the exclusive description. The spirit 



Li 



22 



2411 . There are several popular 



1611 and parton coalescence model Lai inclusively 



of quark combination has been extended to various transport, variation and statistic methods of hadron production in 



relativistic heavy ion collisions 126 



30] 



In this paper, we use SDQCM to treat the initial production of various hadrons. Of all the "on market" combination 
models, SDQCM is unique for its combination rule which guarantees that mesons and baryons exhaust the probability of 
all the fates of the (anti)quarks in deconfined color-neutral system at hadronization. The main idea of the combination rule 
is to line up the (anti)quarks in a one-dimensional order in phase space, e.g., in rapidity, and then let them combine into 



initial hadrons one by one according to this order |]20|. Three (anti)quarks or a quark-antiquark pair in the neighborhood 
form a (anti)baryon or a meson, respectively. The exclusive nature of the model make it convenient to predict the K*° and 
E* production on the basis of the reproduction of the yields and momentum spectra of various stable hadrons. 

In the meson formation, the relative formation probability of the lowest lying vector meson (V) to pseudo-scalar meson 
(P) with the same valance quark composition is tuned by the parameter V/P ratio which can not be given from first 
principles. Spin counting arguments would suggest a 3:1 mixture between vector and pseudoscalar mesons. In Lund 
string fragmentation, the V/P ratio is taken to be 1 for light mesons and 1.5 for strange mesons, which is based on wave 



function overla p arg uments 03 111 . The resulting K* /K yield ratio is about 0.56 which is calculated using PYTHIA 6.4 with 

rT n 

default settings [32J]. The data of K*°/K~ yield ratio in pp collisions at RHIC energies is about 0.35 13]. This will lead to a 
rough estimation of the V/P ratio by the relationship V/(V+P) =0.35 after taking into account K* decay, and the resulting 
V/P ratio is about 0.5. The SDQCM fit of the value of K*°/K~ in pp reactions gives V/P ratio of 0.45. The choice of 
V/P ratio influence directly the predicted abundance of K*° at hadronization, and thus influence the identification of the 
magnitude of rescattering effects in hadronic stage in explaining the experimental data. Here, we give the prediction of 
K*° at different V/P values to test its influence quantitatively. 

The left panel (a) in Figfl]is the pj spectra of various stable hadrons at midrapidity in central Au+Au at y/s^N - 200 



GeV Symbols are the experimental data and lines are results of SDQCM in Ref. l23ll in which the pj spectra of their 
anti -particles are also studied in detail. V/P ratio is taken to be 3 in the calculation. The input of model is the p? spectra 
of light and strange quarks just before hadronization, i.e. fq(pr) and fs(pr), which are fixed by the data of n~ and K~. 
Clearly, the transverse momentum spectra of these stable hadrons can be self-consistently explained by two quark pj 
spectra via combination. This means that constituent quark degrees of freedom play a dominated role in the production 
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FIG. 1: Panel (a): The p T spectra of stable hadrons at midrapidity in central Au+Au collisions at \/s NN = 200 GeV. Panel (b): The p T 
spectrum of A"* . Symbols are the experimental data Ipl l33l - l35n and lines are the results of SDQCM. Panel (c): Yield ratio of K*°/K~ at 
midrapidity. Open squares are results by hadronization for different V/P ratios. Experimental data of pp collisions and central Au+Au 
collisions are shown as dashed line and band area |5|], respectively. Panel (d): Elliptic flow v2 of A"* as the function of p T in minimum 
bias Au+Au collisions at ysiwv = 200 GeV. Symbols are the experimental data |5|] and the solid line is the result of SDQCM. Dashed 
lines marked by n = 3 and 4 are the guidance of the consequence of aggravating regeneration effects in hadronic stage, which is from 
Ref. 
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f. El- 



of these thermal hadrons. 

The middle panel (b) in FigQ]shows midrapidity pj spectrum of K*° produced by hadronization based on the left panel 
results and parameters (quark spectra and V/P ratio). The V/P ratio is sensitive to the K*° yield and is relatively less 
sensitive to the slope of AT* pj spectrum. The K*° spectrum at other V/P values, e.g. 1 and 0.45, are all presented. Here, 
we stress that even at V/P=l and 0.45 the model can still reproduce the pj and rapidity spectra of stable hadrons shown 
left. It is because the V/P ratio does not alter the nature of hadron formation, and it just changes the decay contributions 
of resonance to stable hadrons. We find that the slope of AT* pj spectrum in 1.5-3.0 GeV/c is roughly consistent with the 
data but in the low pj region the directly produced AT* by combination is above the data even for V/P=0.45. 

The right panel (c) in Figfl]is the yield ratio N(K*°)/N(K~) at midrapidity for different V/P values. Because kaon and 
AT* both contain same valence quarks, their ratio roughly cancel the effect of strangeness enhancement in relativistic heavy 
ion collisions and thus is sensitive to the mechanism of hadron production. The data of N(K*°)/N(K~) is 0.2 ±0.03 ±0.03 
in central Au+Au collisions and is 0.34 ± 0.01 ± 0.05 in minimum bias pp reactions at s/srw = 200 GeV 15J, respectively. 
We find that the calculated ratio N(K*°)/N(K~) at V/P =3 and 1 are higher than the data of Au+Au collisions which is 
shown as the band area. The result of V/P=0.45 is consistent with the data of pp reactions shown as dashed line but is 
also higher than the data of Au+Au collisions. The over-prediction of hadronization in the low pj region indicates the 



5 

necessity of hadronic rescattering for the suppression of K*° yield. 

To further manifest the hadronization effect in final state observation, we present in panel (d) in FigQ]the elliptic flow 
v2 of K*° at midrapidity in minimum bias Au+Au collisions at -\/Tnn = 200 GeV. This result is calculated using the 
extracted quark v2 in the previous study of v2 of various stable hadrons in Ref. 13611 . V/P ratio does not influence the 
predicted v2 which is irrespective of particle abundance. One can see that the calculated K*° v2, shown as solid line, is in 
well agreement with the experimental data. The dashed lines are the parameterization of hadronic v2 as the function of 



number of constituent quarks n in Ref. |37f|. n — 3 and 4 are the guidance of the consequence of aggravating regeneration 
effects in hadronic stage. The agreement of our result with the data means that the K*° mesons observed experimentally 
mainly come from the hadronization. 

III. HADRONIC RESCATTERING EFFECTS ON K*° PRODUCTION 

The subsequent hadronic rescattering stage after hadronization is simulated by A Relativistic Transport (ART) model 



113 



18H which includes baryon-baryon, baryon-meson, and meson-meson elastic and inelastic scatterings. We apply the 
code in AMPT event generator V2.25t3 IB8I1 which has provided proper extension to higher RHIC energies from original 
SPS energies. The time span of hadronic stage is important for the hadronic rescattering effects such as the suppression 
magnitude of K* a yield. In transport theory, the rescattering time should be infinite in principle. In practice, one would like 
to choose a finite rescattering time, and this choice, to some extent, is similar to the decouple criteria of kinetic freeze-out 
in hydrodynamic theory of heavy ion collisions. In this work, we treat the rescattering time as an adjustable parameter 
and study how long the time span of hadronic phase is favored by the experimental data when the initial produced K*° is 
fixed. 

Figf2]top panel shows the survived ratio of K* abundance N(K*y /N(K*)""" al after hadronic rescattering stage as 
the function of rescattering time t res . Here N(K*)' m " al is the number of K* just after hadronization and N(K*y lnal is the 
number of K* that can be experimentally reconstructed after hadronic re scatterings. K* mesons are reconstructed from 
their hadronic decay channels using pion-kaon invariant mass analysis, and the survived K* incorporates all effects of 
destruction, signal loss and regeneration. One can see that the number of survived K* almost exponentially decreases with 
the rescattering time. The time of K* number reducing by half is about 20 fm/c. Figf2]bottom panel shows the number of 
survived K*° after hadronic rescattering stage as the function of transverse momentum pj at different rescattering time t res . 
We find that the suppression of K*° caused by hadronic rescattering effects is strong in low pj region. This is qualitatively 



consistent with the result of UrQMD transport model 18 
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FIG. 2: Top panel: The survived ratio of A^'yield as the function of hadronic rescattering time t res in central Au+Au collisions at 
■\fsjw = 200 GeV. Bottom panel: The survived K*° as the function of p T at different t res . 

Taking into account of effects of hadronic rescatterings on the yield suppression and spectra distortion, we obtain the 
yield and pj spectrum of final state K*° comparable to experimental data. The top panel in Figj3]shows the yield of final 
survived K*° at midrapidity as the function of rescattering time t res in central Au+Au collisions at -\/snn = 200 GeV. The 
directly produced K*° for different V/P ratios at hadronization is taken to be the starting point of the hadronic rescattering 
stage. For V/P=3 the rescattering time needed to reproduce the experimental data is about 60 fm/c while for V/P=l the 
needed time is about 30 fm/c. For V/P=0.45 favored by pp N(K*°)/N(K~) data the rescattering time is about 13 fm/c. This 
value is consistent with the typical lifespan of 13 ± 3 fm/c in UrQMD simulation of hadronic rescattering stage |9J. The 
estimation of time span between chemical and thermal freeze-out by a thermal model using an additional pure rescattering 



3911 . The bottom panel in Figj3] shows the pj spectrum 



phase is 2.5^[ 5 fm/c from the analysis of K*°/K~ yield ratio 11121 
of final state K*° at midrapidity at three V/P values with the corresponding rescattering times for yield reproduction in 
central Au+Au collisions at 200 GeV. We find that the strong suppression of low pj K*° in hadronic rescattering process 
offset the over-predication of initial K*° by hadronization shown as the band area, and this leads to an obviously improved 
description of pj spectrum of K*°. 

The K*° data at different collision energies enable the further test of hadronization mechanism and the investigation 
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FIG. 3: Top panel: The final state A"* yield at midrapidity as the function of hadronic rescattering time t res in central Au+Au collisions 
at -\js NN = 200 GeV. The symbols connected with lines are calculation results corresponding three differen V/P ratios at hadronization. 
The experimental data are shown as the band area |5||. Bottom panel: The p T spectra of final state K* for three V/P values with the 
corresponding rescattering times for yield reproduction in central Au+Au collisions at 200 GeV. The experimental data are shown as 
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symbols 



The band area shows the result of K" just after hadronization. 



of the energy dependence of hadronic rescattering effects. In FigJU we present the pr spectra of various stable hadrons 
and K*° resonance at midrapidity in central Au+Au collisions at ysjvjv - 62.4 GeV. Symbols are the experimental data 
and lines are the results of SDQCM. The pj spectra of constituent quarks at hadronization are taken to be the thermal 
exponential pattern exp{-mT/T s ). The slope parameter T s is taken to be 0.31 GeV for strange quarks and 0.29 GeV for 



light quarks, which are smaller than those at 200 GeV IJ23II . The numbers and rapidity spectra of constituent quarks and 



antiquarks are fixed by the experimental data of pion and kaon rapidity spectra 04211 . One can see that, similar to Au+Au 
200 GeV, the model results of various stable hadrons are in well agreement with the data. Then we can predict the pj 
spectrum of initial K*° just after hadronization, and the results are presented in the right panel in Fig0] The degree of the 
agreement between the model result and experimental data is also similar to that in Au+Au 200 GeV. The calculated K*° 
yield densities in low pj region exceed the data and become to close to the data as pj rises to 2-3 GeV/c. This indicates 
the influence of the hadronic rescattering on K*° production is still significant at intermediate RHIC energy. 

The top panel in Fig|5] shows the yield of final survived K*° at midrapidity as the function of rescattering time t res in 
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FIG. 4: The p T spectra of stable hadrons (left panel), their anti-hadrons (middle panel) and A" resonance (right panel) at midrapidity 



in central Au+Au collisions at -\js NN = 62.4 GeV. Symbols are the experimental data 140114111 and lines are results of SDQCM. 
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FIG. 5: Top panel: The final state K* yield at midrapidity as the function of hadronic rescattering time t ms in central Au+Au collisions 
at -\/s NN = 62.4 GeV. The symbols connected with lines are calculation results corresponding three differen V/P ratios at hadronization. 
The experimental data are shown as the band area |5j|. Bottom panel: The p T spectra of final state K* for three V/P values with the 
corresponding rescattering times for yield reproduction in central Au+Au collisions at 62.4 GeV. The experimental data are shown as 
symbols [5J]. The band area shows the result of K* just after hadronization. 



central Au+Au collisions at ysMv = 62.4 GeV. For V/P=3 the rescattering time needed to reproduce the experimental 
data is about 20 fm/c while for V/P=l the needed time is about 14 fm/c and for V/P=0.45 favored by pp N(K*)fN(K~ ) 
data the time is about 5 fm/c. It is found that the needed time at 62.4 GeV is about half of that at 200 GeV. 

The bottom panel in Fig|5] presents the pj spectrum of final state K*° at midrapidity at three V/P values with the 
corresponding rescattering times for yield reproduction in central Au+Au collisions at 62.4 GeV. The strong suppression 
in low p T region in hadronic rescattering process offset the over-predication of initial K*by hadronization shown as the 
band area, and this leads to an improved description of pj spectrum of K*°. 

IV. E* PRODUCTION BY HADRONIZATION AT RHIC 

The E* hyperon has a very short lifetime (4.5 fm/c). The measurement of STAR Collaboration found that the S*/A 
yield ratio at midrapidity in central Au+Au collisions is nearly the same as that in pp collision at -\/snn = 200 GeV 16J]. 
This indicates that the net effects of rescattering loss and rescattering gain in hadronic stage are very small. This poses 
serious constraint on cross sections of various reaction channels of collision gain and collision loss. The small net effect 
of hadronic rescattering has two possibilities. The first is that both the effect of collision gain and that of collision loss 
on X* yield are large but they offset with each other. The pr spectrum of I* might change during hadronic rescattering 
in this case. The second is that two components are all small and pj spectrum of 2* should not be changed dramatically 
compared with the result of hadronization. Here, using the quark spectra fixed by stable hadrons in FigJT] we use SDQCM 
to study the yield and p? spectrum of X* at midrapidity in central Au+Au collisions at y/s^N - 200 GeV. The yield density 
dN/dy of S* + £ at midrapidity just after hadronization is 9.0, and the data of STAR Collaboration is 9.3 ± 1.4 ± 1.2 Q6Q. 
Here S* represents E* + + £*". The result of pj spectrum of E* just after hadronization is presented in left panel in Figj6] 
and is compared with the experimental data. We find that the agreement between hadronization results and the data is 
well. The result favors the second case. We further predict the elliptic flow v2 of S* in minimum bias Au+Au collisions 
at 200 GeV and the result is shown in middle panel in Fig. [6] It is well known that the elliptic flow is sensitive to the 
mechanism of hadron production. The elliptic flow of hadrons formed via quark combination mechanism follows the 
constituent quark number scaling rule. The phenomenological regeneration of S* by A + n — > 2* will result in higher v2. 
To test the effects of hadronic rescattering on Z* production at RHIC energies, we further predict the Z* pj spectrum in 
central Au+Au collisions at 62.4 GeV. The yield density of Z* at midrapidity is 4.3 just after hadronization. 
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FIG. 6: Left panel: The p T spectrum of S* in central Au+Au collisions at ^s NN = 200 GeV. Symbols are the experimental data Ig] and 
the line is the result of SDQCM. Middle panel: The prediction of elliptic flow v2 of E* in minimum bias Au+Au collisions at 200 GeV. 
Right panel: The prediction of p T spectrum of E* in central Au+Au collisions at yjs NN = 62 A GeV. 



V. SUMMARY 



In this paper, we have used SDQCM model for the QGP hadronization and ART model for the subsequent hadronic 
rescattering process to investigate the production of K*° and 2* resonances in central Au+Au collisions at -^snn — 
200 GeV and 62.4 GeV. The initial K*° produced by quark combination mechanism is higher than the experimental 
data in the low pj region and is close to data at 2-3 GeV/c. In the subsequent hadronic rescattering stage, the number 
of K*° that can be reconstructed experimentally exponentially decreases with the increasing rescattering time, and the 
suppression of K*° yield focuses on low pj region, which offsets the over-prediction of hadronization. Therefore the 
quark combination hadronization plus hadronic rescattering can provide a well description of the experimental data of 
K*° production. According to the suppression magnitude of K*° yield, the time span of hadronic rescattering stage is 
estimated. V/P ratio at hadronization is important for the extraction of hadronic rescattering time from the data of K*° 
yield. For V/P =0.45 based on the data of K* a /K~ in pp reaction, the time span of hadronic rescattering stage is about 13 
fm/c at 200 GeV and 5 fm/c at 62.4 GeV. Higher V/P ratio leads to the longer rescattering time. The yield density and pj 
spectrum of E* directly from quark combination hadronization in central Au+Au collisions at 200 GeV is found to be in 
well agreement with the experimental data, which indicates a weak hadronic rescattering effects. To make a further test, 
we predict the v2 of E* in minimum bias Au+Au collisions at 200 GeV and pr spectrum of £* at lower 62.4 GeV for the 
comparison with future STAR data . 
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